mm m 

CELLULAR NEUROSCIENCE 



ORIGINAL RESEARCH ARTICLE 

published: 25 August 2014 
doi: 10. 3389/fncel. 2014. 00246 




Regulation of output spike patterns by phasic inhibition in 
cerebellar granule cells 

Thierry R. Nieus 1t , Lisa Mapelli 23t and Egidio D'Angelo 2 - 3 * 

' Department of Neuroscience Brain Technology, Istituto Italiano di Tecnologia, Genova, Italy 

2 Neurophysiology Unit, Department of Brain and Behavioral Sciences, University of Pavia, Pavia, Italy 

3 Neurophysiology, Brain Connectivity Center, C. Mondino National Neurological Institute, IRCCS, Pavia, Italy 



Edited by: 

Graziella DiCristo, University of 
Montreal, Canada 

Reviewed by: 

Marco Canepari, INSERM, France 
Laurens Bosman, Erasmus MC, 
Netherlands 

'Correspondence: 

Egidio D'Angelo, Neurophysiology 
Unit, Department of Brain and 
Behavioral Sciences, University of 
Pavia, Via Forlanini 6, Pavia, Italy 
e-mail: dangelo@unipv.it 

' First co-authors. 



The complex interplay of multiple molecular mechanisms taking part to synaptic 
integration is hard to disentangle experimentally. Therefore, we developed a biologically 
realistic computational model based on the rich set of data characterizing the cerebellar 
glomerulus microcircuit. A specific issue was to determine the relative role of phasic 
and tonic inhibition in dynamically regulating granule cell firing, which has not been 
clarified yet. The model comprised the excitatory mossy fiber — granule cell and the 
inhibitory Golgi cell — granule cell synapses and accounted for vesicular release processes, 
neurotransmitter diffusion and activation of different receptor subtypes. Phasic inhibition 
was based on stochastic GABA release and spillover causing activation of two major 
classes of postsynaptic receptors, a1 and ot6, while tonic inhibition was based on 
steady regulation of a Cr leakage. The glomerular microcircuit model was validated 
against experimental responses to mossy fiber bursts while metabotropic receptors were 
blocked. Simulations showed that phasic inhibition controlled the number of spikes during 
burst transmission but predicted that it specifically controlled time-related parameters 
(firing initiation and conclusion and first spike precision) when the relative phase of 
excitation and inhibition was changed. In all conditions, the overall impact of a6 was larger 
than that of ot1 subunit-containing receptors. However, a1 receptors controlled granule cell 
responses in a narrow ±10 ms band while a6 receptors showed broader ±50 ms tuning. 
Tonic inhibition biased these effects without changing their nature substantially. These 
simulations imply that phasic inhibitory mechanisms can dynamically regulate output spike 
patterns, as well as calcium influx and NMDA currents, at the mossy fiber — granule cell 
relay of cerebellum without the intervention of tonic inhibition. 
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INTRODUCTION 

Synaptic inhibition controls local microcircuit functions by reg- 
ulating membrane potential and calcium influx in neurons, with 
important consequences on spike generation and synaptic plas- 
ticity (Llinas et al., 2005; Mann and Paulsen, 2007; Lamsa et al., 
2010; Wenner, 2011; Griffen and Maffei, 2014). Since neurons 
can tune their spike patterns on the millisecond scale (Eccles, 
1973; Timmann et al., 1999), it is important to understand how 
the interplay of excitatory and inhibitory synapses regulates spike 
timing. Although recent works have raised attention toward the 
role of tonic inhibitory mechanisms, the role of phasic inhibition 
was less considered. This reflects the difficulty in dissecting the 
different components of inhibition pharmacologically and in acti- 
vating excitatory and inhibitory fibers independently in arbitrary 
phase relationship. 

A prototypical case is posed by the GABAergic inhibitory 
synapses made by Golgi cells on cerebellar granule cells (Eccles 
et al, 1964). Granule cells are normally silent at rest and are 
activated by the mossy fibers. The Golgi cells in turn are acti- 
vated both by mossy fibers and granule cells (Cesana et al., 2013) 



setting up a double feed-back and feed- forward inhibitory circuit. 
The inhibitory process occurs in the cerebellar glomerulus, a spe- 
cialized structure enwrapped into a glial sheet, in which several 
granule cell dendrites are activated by mossy fibers and inhib- 
ited by Golgi cells. The limited diffusion space in the glomerulus 
was shown to favor the establishment of tonic GABA levels and 
various effects of neurotransmitter spillover and cross-talk (Vos 
et al, 1999; Forti et al, 2006; Solinas et al, 2007a,b; Kanichay 
and Silver, 2008). Thus, two inhibitory mechanisms coexist in 
granule cells. Tonic inhibition, which was shown to control the 
granule spike number (Brickley et al., 1996) and transmission 
gain (Mitchell and Silver, 2003). Phasic inhibition, which was 
shown to limit the duration of granule cell responses (D'Angelo 
and De Zeeuw, 2009) and was also implicated in different func- 
tions requiring dynamic network control, including generation 
of granular layer coherent oscillations and resonance (Maex 
and De Schutter, 1998; Dugue et al, 2009; Solinas et al, 2010; 
Gandolfi et al., 2013), induction of mossy fiber - granule cell LTP 
and LTD (Mapelli and D'Angelo, 2007; D'Errico et al., 2009), 
and spatio-temporal reconfiguration of granular layer activity 
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(Mitchell and Silver, 2003; Mapelli et al, 2010a). Phasic and 
tonic inhibition are based on different al- and a6-subunit con- 
taining GABA-A receptors and on their differential localization 
and activation by direct release, spillover and tonic GABA lev- 
els (Nusser et al, 1995; Brickley et al., 1996; Rossi and Hamann, 
1998; Hamann et al., 2002; Rossi et al., 2003; Farrant and Nusser, 
2005; Mapelli et al, 2009, 2014). Although much emphasis has 
been given to the role of tonic inhibition in regulating granu- 
lar layer functions, how these mechanisms dynamically regulate 
granule cell spike patterns has remained largely unexplored. This 
is somehow surprising for the granular layer, a neuronal micro- 
circuit which is thought to provide dynamic regulation of spike 
timing (Eccles, 1973; Fujita, 1982; Timmann et al, 1999). 

By exploiting the advanced knowledge on microstructural 
and molecular properties available, we have modeled the impact 
of phasic and tonic GABAergic inhibition on granule cell out- 
put spike patterns in response to coactivation of excitatory and 
inhibitory synapses by a short mossy fiber burst (Pellerin and 
Lamarre, 1997; Hartmann and Bower, 1998; Chadderton et al., 
2004; Rancz et al, 2007; Arenz et al, 2008) and have predicted 
the impact of single-pulse activation of excitatory and inhibitory 
synapses with variable phase lag and arbitrary combinations of 
excitatory and inhibitory synapses (Pellerin and Lamarre, 1997; 
Hartmann and Bower, 1998; Chadderton et al, 2004; Rancz et al, 
2007; Arenz et al., 2008). Model simulations revealed a critical and 
specific role for phasic inhibitory mechanisms in dynamic regula- 
tion of granule cell firing and calcium currents without the need 
of tonic inhibition. 

METHODS 

EXPERIMENTAL METHODS 

Whole-cell patch-clamp recordings 

Patch-clamp recordings in acute cerebellar slices were performed 
as reported previously (DAngelo et al, 1995; Armano et al., 
2000; Mapelli et al, 2009). Briefly, 17- to 23-day-old Wistar 
rats were anesthetized with halothane (Aldrich, Milwaukee, WI) 
and killed by decapitation. Parasagittal 220 u,m-thick acute slices 
were cut from the cerebellar vermis in cold Krebs solution and 
maintained at 30° C before being transferred to a 2-ml record- 
ing chamber mounted on the stage of an upright microscope. 
The preparations were perfused with Krebs solution (2 ml/min) 
and maintained at 30° C with a Peltier feedback device (TC- 
324B, Warner Instruments, Hamden, CT). Krebs solution for 
slice cutting and recovery contained (in mM): 120 NaCl, 2 
KC1, 1.2 MgS0 4 , 26 NaHC0 3 , 1.2 KH 2 P0 4 , 2 CaCl 2 , and 
11 glucose, and was equilibrated with 95% 02-5% CO2 (pH 
7.4). The patch-clamp pipette solution contained (in mM): 81 
Cs 2 S0 4 , 4 NaCl, 2 MgS0 4 , 0.02 CaCl 2 , 0.1 BAPTA, 15 glu- 
cose, 3 ATP-Mg, 0.1 GTP, and 15 HEPES. For experiments 
with current-clamp recording the internal solution contained 
(in mM): 126 potassium gluconate, 4 NaCl, 5 Hepes, 15 glu- 
cose, 1 MgS0 4 , 0.1 BAPTA-free, 0.05 BAPTA-Ca, 3 Mg 2 +-ATP, 
0.1 Na + -GTP. These solutions maintained resting free [Ca 2+ ] 
at 100 nM and pH was adjusted to 7.2. Patch-clamp pipettes 
filled with these solutions had a resistance of 5-8 MS before 
seal formation. Transient current analysis yielded a membrane 
capacitance (Cm) of 3.9 ± 0.3 pF, membrane resistance (Rm) 



of 2.0 ± 0.2 MS, and series resistance (Rs) of 17.1 ± 0.9 MQ, 
(n = 7). The 3-dB cell plus electrode cut-off frequency was fvc = 
(27rR s C m ) _1 = 2.5 ± 0.1 kHz (n = 7). All drugs were obtained 
from Sigma, except BAPTA tetrapotassium salt (Molecular 
Probes, Eugene, OR); GABA-B receptor blockers CGP35348 and 
CGP55845, mGluRs blockers MCPG, CPPG and gabazine (SR- 
95531, Tocris-Bioscience). All experiments were conducted in 
accordance with international guidelines from the European 
Community Council Directive 86/609/EEC on the ethical use of 
animals. 

Experiments with metabotropic receptors blockers 

GABA-B receptor blockers CGP35348 and CGP55845, and the 
mGluRs blockers CPPG and MCPG were used at 100, 1, 300, 
and 500 u,M, respectively. Perfusion of 10 |xM gabazine further 
blocked the GABA-A receptor mediated inhibition. Incoming 
mossy fibers were stimulated with 4 pulses at 100 Hz, deliv- 
ered every 10 s. Initially, no blockers were added to the exter- 
nal solution, the cells were maintained at the holding potential 
of — 65 mV and the granule cells spikes in response to the stim- 
uli were recorded. The same pattern was repeated 30-40 times 
in control, in presence of metabotropic receptors blockers and 
after perfusion of gabazine (added to the previous blockers, 
Figure IB). 

Tonic current estimates 

The amount of tonic inhibition for each recording was esti- 
mated measuring the holding current change (Vh = — lOmV) 
after gabazine perfusion (on average 21.50 ± 2.12 pA n = 6, 
Figure IF). 

Neurotransmission parameter estimates 

The specific synaptic organization of the cerebellar glomerulus 
enables a precise determination of synaptic parameters. EPSCs 
and IPSCs amplitude, recorded at Vh = —70 and — lOmV 
respectively, was used to estimate the number of excitatory and 
inhibitory synapses activated by the stimulus. The number of 
mossy fibers were estimated by the ratio between peak ampli- 
tude of EPSCs and single-fiber responses (sfEPSCs) recorded 
from granule cells in similar experimental conditions (Sola et al, 
2004). With an average EPSC amplitude = -81.0 ± 12.9pA(« = 
7) and sfEPSC amplitude = -25.8 ± 3.3 pA, the EPSC/sfEPSC 
amplitude ratio yielded an average of 3.14 activated fibers per 
granule cell. The number of activated Golgi cell axons was esti- 
mated by the ratio between peak amplitude of evoked IPSCs 
(elPSC) and spontaneous IPSCs (sIPSC) in the same granule cells. 
The sIPSCs depend on Golgi cell autorhythmic activity and are 
generated at single connections (Hamann et al., 2002; Rossi et al, 
2003; Mapelli et al, 2009). On average, elPSCs amplitude was 
48.5 ± 17.8 pA, while sIPSCs amplitude was 16.3 ± 1.5 pA, so that 
elPSC/sIPSC amplitude ration yielded 2.97 GoC axons activated 
per granule cell (see Mapelli et al., 2009). The parameter space 
explored in the model was based on the estimates of the number 
of input contacts (i.e., 1-4 MFs and 0-4 GoC fibers). The com- 
parison between experimental data and model estimates therefore 
was shown for the 5 granule cells that showed up to 4 GoC and 
MF inputs (cf. Figure IE, see Figures 4, 5 and Table 1). 
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FIGURE 1 | The impact of GABA-A receptors on granule cell excitation. 

(A) Schematic drawing of the cerebellar glomerulus. The granular layer of the 
cerebellar cortex receives input from mossy fibers (MF), making synaptic 
contacts with granule cells (GrC) and Golgi cells (GoC) in specialized 
structures called glomeruli [light gray area). GrC axons (parallel fibers - pf) 
contact GoC dendrites. Therefore, GoCs inhibit GrCs through a feedforward 
(MF-GoC-GrC) and a feedback (MF-GrC-GoC-GrC) pathway [arrows) (modified 
from Mapelli et al., 2014). (B) The traces show patch-clamp recordings of GrC 
activity following a brief burst of MF stimulation (100 Hz, arrows). 
Glutamatergic and GABAergic metabotropic receptor blockers (MCPG, CPPG, 
CGP35348, and CGP55845) modify GrC response. The subsequent blockage 
of ionotropic GABA-A receptors (gabazine) causes an increase in GrC 



excitability and spike generation. (C) The plots show the number of spikes 
(ns), time window (tw), time to first spike [ttfs) and first spike standard 
deviation (fssd) in control icon) and after metabotropic receptors blockage 
(metR) and GABAergic ionotropic receptor blockage (ionR) in 7 GrC 
recordings performed as illustrated in (B). (D) Average EPSC and IPSC 
recorded from the same cell as in (B). Note the different holding potential 
used to isolate the two currents. The gray lines indicate the amplitude of 
unitary EPSCs and unitary IPSCs. (E) The plot reports the estimated number 
of MF and GoC connections active in each one of the 7 experiments covered 
in this figure. (F) The plot shows that the tonic current measured after 
gabazine perfusion in each experiment does not correlate with the number of 
GoC connections (linear fitting slope —0.65, R 2 = 0.137). 



MODELING 

In this paper, a computational model of neurotransmission and 
excitation in the cerebellar glomerulus is developed and a series 
of routines to investigate the parameter space is implemented. 



A single-compartment model of the granule cell [(D'Angelo 
et al., 2001) updated in Nieus et al. (2006)] was imple- 
mented with advanced representations of inhibitory synapses 
and run using the NEURON simulator (NEURON version 7.3; 
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Table 1 | The table summarizes the number of excitatory (n MF) and 
inhibitory (n GoC) synapses activated during each experimental 
recording (exp 1-7). 





n MF 


n GoC 


expl 


3 (3.17) 


1 (0.89) 


exp2 


1 (0.64) 


10 (9.82) 


exp3 


3 (2.63) 


2 (1.65) 


exp4 


4 (3.57) 


1 (0.91) 


exp5 


3 (3.21) 


4 (3.90) 


exp6 


4 (3.75) 


3 (3.02) 


exp7 


5 (5.02) 


1 (1.06) 



The number of excitatory synapses was estimated by comparing the average 
amplitude of eEPSCs with the unitary EPSC amplitude (reported in Sola et al, 
2004), while the number of inhibitory synapses was estimated by comparing the 
average elPSCs amplitude with the sIPSCs amplitude in these same recordings 
(see Methods). The estimated number of active synapses is indicated in brack- 
ets, while the unitary value was used to compare the experimental results to 
the model. 

Hines and Carnevale, 2001). The model was coupled to a record- 
ing electrode to reproduce realistic current-clamp conditions and 
simulations were compared to cellular responses corrected for 
10 mV liquid-junction potential. The ionic mechanism were the 
same as reported in Nieus et al. (2006) and Solinas et al. (2010). 
The granule cells are extremely compact also when channels are 
open in the dendrites during synaptic transmission. Calculations 
based on cable equation (D'Angelo et al., 1993, 1995; Rossi et al, 
1994) and tests using bio-realistic compartmental models (Silver 
et al., 1992; Diwakar et al., 2009) showed that dendritic electro- 
tonic length was around L — 0.04 and the maximum current and 
membrane potential loss was <5%. Therefore, the results would 
not be different using single or multi-compartment models. 

Major aspects of the inhibitory mechanisms used as the base for 
modeling 

In the cerebellar glomerulus, GABA released from the Golgi cell 
presynaptic terminal can act directly on receptors facing the 
postsynaptic site (fast direct component) or indirectly on extrasy- 
naptic receptors through spillover from neighboring sites (slow 
indirect component) (Rossi and Hamann, 1998; Hamann et al., 
2002; Rossi et al, 2003; Mapelli et al, 2009). Direct release from 
the presynaptic terminal activates al GABA- A receptors and, in 
part, a6 GABA-A receptors. While al receptors are located in 
the PSD, a6 receptors are located preferentially at extrasynap- 
tic sites (but in part also in the PSD, Nusser et al, 1995). Fast 
kinetics of the direct component are explained by rapid gating 
of al -containing GABA-A receptors. Conversely, slow kinetics of 
the indirect component are explained by slow gating of a6 recep- 
tors largely due to spillover from neighboring contacts (Rossi and 
Hamann, 1998; Hamann et al., 2002; Rossi et al, 2003). In addi- 
tion, ambient GABA at sub-micromolar concentration activates 
high affinity GABA-A receptors generating tonic inhibition. 

Mathematical model of inhibitory synaptic transmission 

A mathematical model of Golgi cell — granule cell neurotrans- 
mission was developed considering available knowledge on the 



synapse and on GABA-A receptor kinetics (Cherubini and Conti, 
2001; Farrant and Nusser, 2005). The Golgi cell — granule cell 
IPSCs are generated by different GABA-A receptor subtypes: ai 
subunit-containing receptors contribute mainly to early IPSC 
activation and determine the IPSC peak, while subunit- 
containing receptors sustain the IPSC along its decay phase (Tia 
et al, 1996; Rossi and Hamann, 1998; Brickley et al., 2001). 
Differential receptor activation reflects a range of GABA affinities 
for the different subtypes. In particular, ai-subunit containingre- 
ceptors have EC50 in the 10-100 |xM range, ag|32/3Y2 receptors 
have EC50 in the |xM range, and a^2/3^> receptors have EC50 
in the 10 nM range (Tia et al, 1996; Nusser et al., 1998; Hadley 
and Amin, 2007). Differential receptor activation also reflects 
their subcellular localization. While ai receptors are localized in 
the synaptic cleft, 06^2/38 receptors are hundreds of nanometers 
away from the postsynaptic densities (Nusser et al., 1998), where 
they are most likely activated by low ambient GABA concentra- 
tions and by GABA spillover (10 nM to few |xM, (Farrant and 
Nusser, 2005)). The agf52/3Y2 receptors are localized in the synap- 
tic junction (Nusser et al., 1998). The simulation of GABAergic 
neurotransmission in the glomerulus required explicit represen- 
tations of GABA release and diffusion as well as of the kinetic 
properties of different GABA receptor subtypes. 

Mechanisms of elPSC generation. GABA release was imple- 
mented through a phenomenological 3 -state model (ref. 
Markram and Tsodyks) that well reproduces synaptic dynamics 
on short time scale, as already demonstrated at the mossy-fiber 
to granule cell synapse (Nieus et al., 2006). The scheme was gov- 
erned by recovery rates and release probability from the reserve 
pool and its output was GABA release. The neurotransmitter 
concentration, [GABA], was modeled by summing a pulse to a 
diffusive term: 

[GABA] = [GABA]p + [GABA] D (1) 

[GABAjp activates the receptors facing the releasing site, while 
[GABA]d accounts for GABA spillover in the glomerulus (Rossi 
and Hamann, 1998; Nieus et al., 2006). Neurotransmitter concen- 
tration at a distance r from the releasing site was described by the 
diffusion equation: 

M _=ii_ 

[GABA] D = - e 4 ' p eff (2) 
4rmD eff 

where M represents the amount of released neurotransmitter 
molecules, D e ff is the effective diffusion coefficient of GABA 
molecules, h is the cleft thickness and r is the distance from 
which the neurotransmitter diffuses. Given the strong similar- 
ity of GABA and glutamate molecules (which differ only for 
a carbossilic group), the diffusion coefficient of GABA was set 
equal to that of glutamate (D e ff = 0.221 ixm 2 /ms; Nielsen et al., 
2004; Nieus et al, 2006). From morphometry, h was set at 50 nm 
(Hamori and Somogyi, 1983; lakab and Hamori, 1988). The 
parameters M and r in Equation 2 could not be estimated inde- 
pendently and were obtained from fitting (see below and Table 2). 
To account for the coexistence of direct and indirect components 
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Table 2 | Model parameters obtained from elPSC train fittings with 
the model. 



Control CGP55845 



G$ ax (pS) 1323.25 ± 310.99 

G^ 6 ax (pS) 491.58 ± 112.32 

/\(n°molec) 10573 ± 3324 

r(n,m) 1.07 ±0.54 

Tflfc(ms) 38.7 ± 38.6 

z FAC {ms) 0 

p 0.42 ± 0.08 0.67 ±0.09 



Q*max and G"max are tr, e maximum conductance of a; andae receptors, A is the 
number of diffusing molecules and r is the diffusion distance, t rec and z FAC 
are the time constants for synaptic vesicle recovery and facilitation, p is release 
probability. Note that since in all fits the facilitation time constant remained pretty 
low, tfac w as fixed to 0 without affecting the quality of the fit itself. The data, 
except for p, are in common since control and CGP55845 traces were fitted 
simultaneously and p was the only independent parameter. Data are reported as 
mean ± s. e.m. (n= 5 cells). The model elPSC train fittings and states are shown 
in Figures 2, 3. 

we implemented two separate GABA-A receptor kinetic schemes. 
The ai- and c*6- receptor responses were reproduced using the 
same structural kinetic scheme (see Figure 2A) while the kinetic 
rates of the two receptor subtypes were determined by optimiza- 
tion to data taken from the literature (Figures 2 A,B). The IPSC 
was then reproduced as: 

IPSC = elPSCai + eIPSCa 6 (3) 

The ot6 component mostly sustains the elPSC tail after the 
ai response has almost completely decayed (cf. Wall, 2002). 
Diffusion also contributes to protract the elPSC through the ai 
component. The GABA channels were given a reversal poten- 
tial Vci = —65 mV and the tonic component was modeled as a 
background conductance (D'Angelo et al, 2001). 

GABA-A receptor kinetics. The GABA-A receptor models are 
characterized by two desensitized, three closed and two open 
states (Barberis et al., 2000; Jones et al, 2001). These models 
are based on IPSC reconstruction from out-side-out patch-clamp 
experiments. In whole cell in slices, IPSCs show faster decay than 
in patches, probably reflecting an increased desensitization caused 
by phosphorilation of the GABA receptor (Brickley et al., 2001). 
Interestingly, parallel investigations (Petrini et al., 2011) have 
demonstrated that the discrepancy between slice and out-side-out 
experiments might also be determined by the neurotransmitter 
diffusion in the synaptic cleft that can also shape the IPSC time 
course. 

Recently, a multistate model was proposed that accounts for 
out-side-out and slice IPSC responses (Pugh and Raman, 2005). 
The Pugh-Raman model suggests that increasing the transition 
rate into the fast desensitizing state ("DA2f" in the scheme of 
Figure 2A) can convert out-side-out into slice channel behavior. 
The ai receptor model was obtained by fitting the Pugh-Raman 
model to sIPSCs generated by the GABA pulse. The receptor 
model was obtained by fitting the Pugh-Raman model to slow 



elPSCs activated by GABA spillover (see Figure 2 in Tia et al., 
1996; Rossi and Hamann, 1998; Wall, 2002). Fittings yielded an 
EC50 ~ 1.5 (iM placing the a6-model in the activation range for 
016P2/3Y2 receptors. Then, ai and ct6 receptor currents combined 
to generate the composite elPSC. The kinetic rate constants of ai 
and 0(6 receptors are reported in Table 3. 

Simulation methods. The model was written in NEURON based 
on a previous reconstruction of the granule cell (for details see 
D'Angelo et al., 2001; Nieus et al., 2006) and is available on http:// 
senselab.med.yale.edu/modeldb/default.asp. The model allowed 
to determine neurotransmission parameters through a fitting 
procedure carried out in two steps. First, single elPSCs were fit- 
ted to determine GABA diffusion parameters (except for h and 
D e ff, which were reported from literature). Then, the whole elPSC 
train was fitted on data from our previous experimental investiga- 
tion (Mapelli et al, 2009) yielding parameters related to synaptic 
dynamics including the vesicle recycling and facilitation time 
constants (trec and tfac) and release probability (p). 

Stochastic model. In order to test the GRC responsiveness 
in a more natural context we introduced synaptic stochastic 
mechanisms at the excitatory and inhibitory synapses. We 
modeled stochastic pre-synaptic release of the excitatory and 
inhibitory phasic components (i.e., the ai receptor). The stochas- 
tic quantal release model was implemented as in Arleo et al. 
(2010) by splitting each synapse into a fixed number of indepen- 
dent releasing sites (RS) and independent post synaptic densities 
(PSD). The presynaptic RSs and the postsynaptic PSDs were 
endowed with the same kinetic scheme as in the deterministic 
model. The synaptic conductances were rescaled by the num- 
ber of RSs. Each stimulus, e.g., presynaptic AP, is broadcast to 
the RSs and an uniform [0,1] random number (rnd) is gener- 
ated at each RS where a fixed concentration of neurotransmitter 
(Tmax) is released when the corresponding rnd falls below the 
actual releasable resources (Y). This stochastic model yields the 
deterministic response when averaged over multiple trials (not 
shown). 

Model of the glomerular microcircuit and simulations 

In order to test the responsiveness of the granule cell (GC) in a 
realistic context, a GC microcircuit was built that included the 
model of the GC (D'Angelo et al, 2001) as well as its excitatory 
(Nieus et al., 2006) and inhibitory synapses (developed in this 
manuscript). At excitatory and inhibitory synapses three indepen- 
dent releasing sites were included (Sola et al., 2004; Mapelli et al, 
2009; Arleo et al., 2010), each one comprising the same stochas- 
tic release mechanisms. The GC received a maximum of four 
excitatory mossy fibers (MF) and four inhibitory Golgi (GoC) 
axons. These average numbers were adopted to limit the param- 
eters space and the case of higher numbers of synapses (e.g., 10 
inhibitory synapses in Figure IE) or releasing sites (e.g., see Sola 
et al., 2004; Saviane and Silver, 2007; Mapelli et al, 2009) were not 
accounted for. 

Two stimulation protocols were used to simulate the impact 
of the different inhibitory mechanisms on the granule cell 
output. First, a high frequency input (100 Hz, 4 pulses) was 
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FIGURE 2 | Modeling GABA receptor-mediated responses in the 
cerebellar glomerulus. (A) Schematic representation of the kinetic scheme 
for m and ae subunit-containing receptors (modified from Pugh and Raman, 
2005). The traces [continuous lines) show responses to brief and long 
1 mM GABA pulses (taken from Tia et al., 1996). The simultaneous fitting 
of the decay phase of the ag response to a fast (1 ms) GABA 
neurotransmitter pulse (1 mM) and of the almost constant response to a 
long (200 ms) GABA neurotransmitter pulse (1 mM) is shown {dashed lines) 
(see Methods and Table 3 for details). (B) The simulated elPSC generated 
by the model at OmV in response to a 1 ms GABA pulse shows amplitude 
and kinetics reproducing experimental elPSCs. The contribution of the on 



and <»6 components in the elPSC is shown, revealing that the at receptors 
account for the direct response while the ci6 receptors account for the 
indirect response [dashed lines). (C) The simulations show the impact of 
changing pi (from 0.42 to 0.64) on the GABAergic response to a train of 
impulses at 100 Hz. The pi increase changes the profile of GABA 
concentration (in particular, the first peak of GABA concentration is 
enhanced) modifying the kinetic states of ai and a6 responses (O, C, and 
D indicate the sum of all the open, closed and desensitized states for each 
receptor type). The ai is more affected than as response by pi, since as 
receptors reacts slowly and tends to saturate. The associated changes in 
the ai and <J6 currents are shown at the bottom. 



used to activate the excitatory synapses, while the inhibitory 
synapses were activated by the same input train delayed by 
5 ms to mimic synaptic and integration delays in the feed- 
forward inhibitory MF-GoC-GC pathway. The results obtained 
with the first protocol are reported in Figures 4-6. Secondly, sin- 
gle shock stimuli were delivered to the excitatory and inhibitory 
synapses with a variable phase-lag (from —100 to +100 ms). 
The results obtained with the second protocol are reported in 
Figures 7, 8. 

To the ease of comparison with experimental results we mon- 
itored set of parameters such as the number of output spikes 
(sn), the timing of the first spike (ff/s), the first spike standard 
deviation (fssd, defined as the standard deviation of the mea- 
sure of ttfs) and the time window of the output spikes (tw, 
defined as the time difference between the last spike and the 
first stimulation). Reliable estimates of the measured parame- 
ters were then obtained by running 500 runs for each tested 
condition. 



Table 3 | Kinetic rate constants used to simulate a\ and ae GABA-A 
subunit-containing receptors in cerebellar granule cells. 





«1 


«6 


k on (mM -1 -ms -1 ), k off (ms- 1 ) 


20, 6 


54.8, 0.31 


ai (ms -1 ), bi (ms -1 ) 


0.06, 0.03 


0.06, 0.03 


a2(ms~ 1 ), b2(ms _1 ) 


0.4, 10 


0.4, 10 


di (ms -1 ), n (ms -1 ) 


3.3e-4, 7e-4 


0.86, 0.04 


d2(ms~ 1 ), r2(ms _1 ) 


1.2, 6e-3 


2.7, 0.43 


d3(ms _1 ), r 3 (ms -1 ) 


15, 3.75 


15, 7.41 


di d2(mM -1 ■ ms -1 ), n r2(ms -1 ) 


15, 0.007 


24.2, 0.09 



The corresponding simulated traces and kinetic scheme are shown in Figure 2 



Data analysis 

All data were analyzed with custom codes developed in 
Python (www.python.com) based on the Scipy (www.scipy.com), 
Numpy (www.numpy.com) and MatPlotLib (www.matplotlib. 



Frontiers in Cellular Neuroscience 



www.frontiersin.org 



August 2014 | Volume 8 | Article 246 | 6 



Nieus et al. 



Phasic inhibition in cerebellar granule cells 



com) modules. The analysis allowed to quantify the number of 
spikes elicited, the time window, the time of the first spike and its 
variability (i.e., first spike precision) in the tested E/I conditions. 
In a set of simulations we also recorded the Calcium and NMDA 
current traces and computed the mean values over a 130 ms long 
time window. 

A compact representation of the E/I space (e.g., Figure 4B) 
was obtained by reporting the excitatory and inhibitory strengths 
on Cartesian coordinates using a color-code map, in which 
a color scale corresponded to the intensity of the measured 
parameter. The excitatory, inhibitory strengths (Se, Si) were 
quantified as the sum of the active fibers (MF/GoC) and their 
release probability (Pel Pi), such that Se = MF + Pe and 
Si = GoC + Pi. 

The Se/i was calculated on sub-regions of the E/I space by 
averaging over the active fibers involved and the corresponding 
release probabilities (e.g., Figure 4B). Each sub-region was named 
after the number of active fibers involved: for example, £/7(l/l) 
corresponds to 1 excitatory and 1 inhibitory fiber. Se/i yields 
a comprehensive quantification of the changes of the measured 
quantities at variable E/I levels. 

In order to describe the shape of Se/i data distributions at 
variable phase-lag between synaptic excitation and inhibition, 
the data were fitted with a Lorentz equation y(6) of the form 
(Gandolfi et al., 2013): 



2- A 



x 4 ■ (6 - 6 max ) 2 + co 2 



(4) 



where, A is the underlying area, 9 max is the peak frequency, and 
co is the width at y{0 ma %)l2. Lorentzian fittings allowed to find the 
resonance phases (RF = 0 max ) and to calculate RA = y(9 max )/co, 
which corresponds to the quality factor Q in resonance literature 
and characterizes a resonator's bandwidth relative to its center 
delay. The goodness of fit was assessed by calculating the squared 
correlation factor, R 2 . Q is normally used in signal processing the- 
ory and describes the tendency of a signal to show resonance vs. 
monotonic roll-off. 

The E/I response space representation also distinguishes 
between different settings of the GABA tonic current when the 
phasic component is zero. 

Entropy calculation 

To estimate the redundancy of the parameters (sn, tw, ttfs) 
on the E/I space we calculated the coding fraction (Borst and 
Theunissen, 1999; Sadeghi et al., 2007) defined as: 



CF P 



MI P 
HiS) 



(5) 



where H(S) is the entropy of the E/I input space and MI quantifies 
the mutual information of the measured parameters P (P = sn, 
tw, ttfs). Each Se/i configuration of the E/I space was regarded as 
an independent stimulus. The cardinality of the stimuli is there- 
fore given by: 4 (MF) x 9 (Pe) x 5 (GoC) x 9 (Pi) = 1620. 
CFp measures the fraction of Se/i that can be distinguished based 
on the observation of the whole output generated by the model 
(i.e., 500 trials for each configuration) and is normalized between 



0 (no separation, 100% redundant) and 1 (perfect separation, 
0% redundant). CF calculation was performed with the Python 
module "pyentropy" (Ince et al., 2009) that provides appropriate 
procedures to robustly estimate MI. 

RESULTS 

In this paper we have developed a computational model of the 
cerebellar glomerular microcircuit (Figure 1A) and we have sim- 
ulated its functionality exploring the response space governed 
by synaptic inhibition. Simulation results were compared with 
experimental responses to mossy fiber burst stimulation. 

THE EFFECT OF INHIBITORY NEUROTRANSMISSION ON GRANULE CELL 
RESPONSES TO MOSSY FIBER BURSTS 

Short high-frequency bursts represent a fundamental pattern of 
mossy fiber activity (Chadderton et al., 2004). Bursts co-activate 
granule cells and Golgi cells (D'Angelo et al., 1995; Mapelli et al, 
2014). While granule cells depolarization begins, the progressive 
increase of activity in the inhibitory circuit tends to modulate 
their response. An early component of inhibition passes through 
the mossy fiber — Golgi cell — granule cell circuit (feed-forward 
inhibition), while a late component passes through the mossy 
fiber — granule cell — Golgi cell — granule cell circuit (feed-back 
inhibition) (Cesana et al, 2013; D'Angelo et al., 2013). 

The impact of inhibitory neurotransmission was investi- 
gated by measuring how it modulated granule cell excitation 
following mossy fiber burst stimulation in acute cerebellar slices 
(a condition in which feed-forward inhibition prevails). Since 
metabotropic mechanisms can influence granule cell responses 
(see (Mapelli et al., 2014) for review), ionotropic effects were 
isolated using a cocktail of antagonists for mGluRs and GABA- 
B receptors (CGP35348, CGP55845, MCPG, CPPG) (Figure IB). 
Following application of these blockers, on average the number of 
spikes increased (sn change 39.8 ± 10.5%, p < 0.034, n = 7), the 
discharge time-window increased (tw change 93.4 ± 43.1% p < 
0.034, n = 7) while the time-to-first-spike and first spike standard 
deviation did not vary significantly (ttfs change —2.3 ± 17.5%, 
n = 7, p = 0.79; fssd change -1.17 ± 1.34%, n = 7,p = 0.97) 
(Figure 1C). These results show that metabotropic receptors 
could modify burst transmission to granule cells indicating that 
the data to be compared to our cerebellar microcircuit model sim- 
ulations had to be those obtained in the presence of metabotropic 
receptor blockers. 

Inhibitory neurotransmission between Golgi cells and granule 
cells is mediated by GABA-A receptors, which can be selectively 
blocked by the GABA-A receptor antagonist, gabazine (Rossi 
et al., 2003). This receptor antagonist is normally used in a vari- 
ety of experimental conditions and blocks the phasic and tonic 
inhibitory processes altogether (Hamann et al, 2002; Nusser and 
Mody, 2002; Mitchell and Silver, 2003; Rossi et al, 2003; Mapelli 
et al., 2010a,b; Gandolfi et al, 2013). Following metabotropic 
receptor blockage (see above), the application of 10 \iM gabazine 
significantly increased number of spikes (68.9 ± 41.1%, p < 
0.05, n = 7) and time window (112.6 ± 45.0%, p < 0.04, n = 7) 
but did not significantly change the time to first spike (2.74 ± 
10.83%,p < 0.79, n = 7) or first spike standard deviation (23.4 ± 
80.0%, p < 0.32, n = 7) (Figure 1C). Thus, as a whole, GABA-A 
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receptor-mediated inhibition reduced the intensity and duration 
of granule cell responses without significantly changing the aver- 
age time of occurrence and precision of the first spike. This was 
expected considering that Golgi cell inhibition through the feed- 
forward circuit usually occurs after the granule cell has emitted 
the first spike. 

In each of the recorded granule cells, measurements of EPSCs 
and IPSCs were obtained at the reversal potential of IPSCs and 
EPSCs, respectively (Figure ID). Moreover, the tonic GABA cur- 
rent was estimated by measuring the leakage current before and 
after gabazine application (not shown). We did not found any sig- 
nificant correlations between the tonic current and the estimated 
number of active Golgi cell axons (Figure IF), suggesting that the 
tonic current was regulated by ambient GABA rather than by pha- 
sic neurotransmitter release (for review see Mapelli et al., 2014). 
Specific simulations showed that changes in the tonic current in 
the range observed experimentally did not significantly affect the 
model outcome (<5% difference for spike number, first-spike 
delay, time-window; data not shown). Therefore, in the model 
tonic inhibition was set as a fixed value, not correlated with the 
number of active synapses. 

The number of active synapses was estimated by comparing 
EPSCs and IPSCs with the average single-site responses measured 
previously in identical experimental conditions (Figure ID) (Sola 
et al., 2004; Mapelli et al., 2009). The estimated number of active 
mossy fibers ranged from 1 to 5 and that of active Golgi cell 
synapses ranged from 1 to 10 (Figure IE). Therefore, for each 
recorded granule cell, the most relevant parameters required to 
determine its inhibitory control were available and this allowed a 
direct comparison with simulated data (see Figure 6). 

MODELING INHIBITORY NEUROTRANSMISSION 

A model of cerebellar glomerular inhibition should be able to pre- 
dict the different kinds of synaptic inhibitory responses observed 
experimentally in granule cells. 

mlPSCs 

The minimal GABAergic responses are the miniature synap- 
tic currents, or mlPSCs, which correspond to release of single 
neurotransmitter quanta (Mapelli et al., 2009). 

sIPSCs 

The spontaneous multiquantal GABAergic responses, or sIPSCs, 
are caused by spontaneous activity in Golgi cells (Dieudonne, 
1998; Rossi and Hamann, 1998; Forti et al, 2006; Solinas et al, 
2007a,b; Mapelli et al, 2009; Brandalise et al, 2012; Hull and 
Regehr, 2012). Since sIPSCs reflect activation at single Golgi con- 
nections, sIPSCs are examples of direct release without spillover 
from neighboring contacts. 

elPSCs 

The evoked IPSCs elicited by electrical stimulation of the Golgi 
cell axon, or elPSCs (Rossi and Hamann, 1998; Mapelli et al, 
2009), are composed by both a direct component due to release of 
GABA from the presynaptic site facing the postsynaptic density 
and an indirect component due to spillover of neurotransmitter 
from neighboring presynaptic sites. The decay of elPSCs is indeed 
slower than that of sIPSCs. 



Inhibitory synaptic transmission in the cerebellar glomeru- 
lus was reconstructed considering the following salient aspects 
(see Methods for details): (i) the vesicular release cycle contained 
an explicit representation of the number of releasing sites and 
of release probability, (ii) the mechanism of vesicle fusion was 
stochastic, (iii) the postsynaptic membrane expressed al and a6 
receptors, (iv) GABA diffusion was explicitly modeled, and (v) 
tonic GABA receptor activation was accounted for. 

The granule cell GABA receptor models were characterized 
by two desensitized, three closed and two open states (Barberis 
et al., 2000; lones et al., 2001) whose transitions were differen- 
tially controlled by GABA transients and slow diffusion waves 
(Petrini et al., 2011). The GABA receptor responses were mod- 
eled using a modified Pugh and Raman (2005) 8-state kinetic 
scheme using the parameters optimized to match responses gen- 
erated by GABA pulses (ai -type; Brickley et al, 2001 ) and GABA 
spillover (ct6-type; Tia et al, 1996; Rossi and Hamann, 1998; 
Wall, 2002) (see Methods for details; Figure 2A; Tables 2, 3). The 
elPSCs were reconstructed as the sum of ai and receptor- 
mediated components (Figure 2B). The time course of oti and 
oi6 receptor-mediated currents and their combinations were com- 
patible with the time course of sIPSCs and elPSCs (Mapelli et al., 
2009). 

The response of the model during repetitive synaptic transmis- 
sion is shown in Figure 2C. The elPSC peak showed depression 
along the train, while the elPSC tail showed the marked temporal 
summation typical of experimental recordings. Clearly, subse- 
quent IPSC peaks were driven by oscillations in the activation 
state of ai receptors, while temporal summation was mostly 
determined by integration of GABA signals through receptors. 
The oi6 receptors activated slowly and, due to their high affinity, 
tended to saturate in a few impulses. Interestingly, changing 
release probability had a remarkable impact on ai receptor- 
mediated peaks but much less so on the receptor-mediated 
current plateau, as further considered below. 

MODELING THE EFFECT OF RELEASE PROBABILITY CHANGES 

A raise in release probability was proposed to explain the changes 
in synaptic responses caused by GABA-B receptor blockage by 
CGP55845 (Mapelli et al, 2009). Here, these changes were sim- 
ulated by increasing the release probability from pi = 0.42 to 
pi = 0.67 (Figure 2C). This transiently enhanced GABA concen- 
tration in the synaptic cleft, but within 4 impulses the concen- 
tration settled toward similar levels with both p values. While 
the ai receptors followed the GABA concentration changes, the 
oi6 receptors developed similar kinetics at both release probabili- 
ties. The slow reaction of ae receptors and their early saturation 
explained the invariance of temporal summation and total charge 
transfer despite the clear changes caused by pi on the ai receptor- 
mediated response. 

In order to test whether a p change could indeed provide a 
global interpretation of IPSC changes during trains, control and 
CGP55845 elPSC trains reported by Mapelli et al. (2009) were 
simultaneously fitted with the model and pi was left free to vary 
(Figure 3A). The fittings yielded pi = 0.42 ± 0.12 in control and 
pi = 0.67 ± 0.13 in CGP55845 (« = 5), in close agreement with 
experimental determinations (n = 5;p = 0.14, paired r-test). The 
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FIGURE 3 | Fitting elPSC trains with the model. (A) Fittings of 
experimental traces in control and after CGP55845 perfusion (taken from 
Mapelli et al., 2009). Both traces are fitted simultaneously with the model 
maintaining the same a-\ and ag receptor density and kinetic schemes, the 
same diffusion parameters and the same presynaptic (tr, xp) parameters. 
Only pi is left free to vary independently in the two traces. The fittings 
properly follow the experimental traces by adapting pi as the only 
independent parameter. (B) The plot shows normalized peak amplitude 
(mean ± MSE) obtained from fittings (n = 5) of elPSC trains (n = 5). Similar 
to real data measurements (see Mapelli et al., 2009), the time course of 
changes after increasing pi demonstrates the increase in the first elPSC 
and the tendency to attain control amplitude in subsequent responses 
(mean ± MSE). The histogram represents charge transfer during the 
simulated elPSC trains, showing no significant difference between control 
and after increasing pi {white bar and black bar; p = 0.55 paired f-test). 



normalized elPSC amplitude showed a significant difference on 
the first pulse but then recovered to control values and the total 
electrical charge transferred during the IPSC train was in good 
agreement with the experimental one and did not change sig- 
nificantly by raising pi (from 1.19 ± 0.20 to 1.46 ± 0.24 pC; 
n = 5; p = 0.31 paired t-test) (Figure 3B). These simulations 
indicate that a raise in release probability can indeed explain a 
consistent set of changes in composite ai-a^ receptors-mediated 
synaptic currents during repetitive neurotransmission suggesting 
appropriate setting of synaptic transmission parameters. 

MODELING THE EFFECT OF INHIBITION ON GRANULE CELL FIRING 

In order to quantify granule cell responsiveness in different acti- 
vation regimens, several series of simulations were run using the 
dendritic microcircuit model (Figure 4A) and the spiking parame- 
ters (number of spikes, time-window, time to first spike and first 
spike standard deviation) were calculated from raster plots and 
PSTHs of the responses (see Methods). An exhaustive exploration 
over all possible excitatory/inhibitory (E/I) synapse combinations 
and release probabilities was performed generating multipara- 
metric color plots (Figure 4B). 

Some regularities emerged from this representation. First of 
all, the parameters varied in a predictable manner with respect 



to the E/I balance: number of spikes and time window appar- 
ently increased by several times while time to first spike and first 
spike standard deviation showed modest changes by raising E/I. 
On a macro-scale, in all graphs an ideal bisectrix separated a low 
response region (top left, low E/I) from a high response region 
(bottom right, high E/I). Moreover, a micro-pattern emerged 
inside each sub-quadrant determined at each E/I by varying 
excitatory and inhibitory release probabilities, pElpI. Thus, in 
principle, by combining different numbers of fibers with differ- 
ent release probabilities, a broad set of response combinations 
could be obtained generating an extended response space. A 
quantitative insight on its properties was obtained by consider- 
ing representative areas. For example, by comparing parameters 
in the quadrants El,2/I3,4 (low E/I) and E3,4/I1,2 (high E/I) 
there was a nearly 10-fold increase in number of spikes, a 3-fold 
increase in time window and a modest changes in time to first 
spike and first spike standard deviation (Figure 5A). 

The inspection of the graphs reported in Figure 4B showed 
that, after switch-off of the phasic components of inhibitory 
synaptic transmission, number of spikes and time window 
increased in all quadrants, while time to first spike and first 
spike standard deviation increased only at low E/I. Interestingly, 
the effect of ot6-receptor switch-off dominated over that of al- 
receptor switch-off in all cases. It should also be noted that the 
effects of changing E/I and pElpI, as well as the effects of al and 
ct6-receptors switch-off, tuned first-spike timing and precision 
over the sub-millisecond time-scale and that the spike num- 
ber varied in discrete units of a few spikes, making the system 
highly sensitive to activity variations. The complete switch-off 
of all forms (phasic and tonic) of GABAergic inhibition (simu- 
lating a gabazine application see Figure 1) caused a proportion- 
ate increase only in spike number and time-window duration 
without changing time to first spike and first spike standard 
deviation. 

Finally, the model showed that raising pi caused a modest 
decrease (<20%) of number of spikes, time window, time to first 
spike and first spike standard deviation. Conversely, raising pE 
remarkably increased number of spikes (<80%) and decreased 
time to first spike and first spike standard deviation (<80%) with 
a limited effect on time window (<20%) (Figure 5B). Therefore, 
the effect of synaptic inhibition on burst parameters were more 
related to the number of active inhibitory synapses and to E/I than 
to the inhibitory release probability. 

EXPERIMENTAL DATA IN THE EXCITATORY/INHIBITORY SPACE 

The simulated parameters were compared to those measured 
experimentally (cf. Figure 1). The points were reported in their 
sub-quadrant following the indications derived from Figure IE 
and located at the corresponding pElpI intersection (Figure 6A). 
In this comparison we used only the 5 granule cells with E, I = 4 
in order to fit into the simulated parameter space. For these 
cells, the pE and pi values closely corresponded to those reported 
previously (Sola et al., 2004; Mapelli et al, 2009) (Figure 6B). 
The parameter variation caused by gabazine was then compared 
to that caused by switch-off of both tonic and phasic synaptic 
inhibition in the model (Figure 6C). The data measured for the 
different parameters showed a variation similar to that of the 



Frontiers in Cellular Neuroscience 



www.frontiersin.org 



August 2014 | Volume 8 | Article 246 | 9 



Nieus et a 



Phasic inhibition in cerebellar granule cells 




FIGURE 4 | Effect of the phasic and tonic inhibition on granule cell firing 
in response to mossy fiber bursts. (A) Exemplar voltage traces simulated 
using MF = 2. Pe = 0.3, GoC = 2, Pi = 0.2 [cf. coordinates (2.3, 2.2) in the 
E/l space in B]. The traces show responses in control and after switch-off of 
ai and ci6 subunit containing receptors. (B) E/l response space for number of 
spikes (sn), time window (tw), time to first spike (ttfs), first spike standard 
deviation {fssd) in control (con) and after switch-off a-i and «6 subunit 



containing receptors (a1-OFF and a6-OFF) and of the GABA tonic current 
(tonic OFF). The E/l space is generated by reporting the number of active 
excitatory (1 < E < 4) and inhibitory (0 < / < 4) synapses on the two axes 
and by using colors to represent the parameter values. For each E/l 
combination, a subquadrant identifies a pE/pl subspace (0 < pE < 1 on the 
x-axis, 0 < pi < 1 on the y-axis). Note that ag subunit-containing receptors 
are most effective in modifying the parameters. 



model (linear fitting with slope 0.95, R 2 = 0.914). It should be 
noted that this data distribution was not significantly different 
from the unitary diagonal (X 2 p < 2.6 e~ 10 ). Finally, it should 
be noted that time to first spike and first spike standard devia- 
tion showed modest changes compared to number of spikes and 
time window, in agreement with experimental results matched 
for the E/I balance (experimental E/I = 2.2 compared to simu- 
lated E/I = 2.6). Therefore, the model was able to capture the 
main quantitative properties of glomerular transmission and the 
impact of synaptic inhibition on granule cell firing. 



between mutual information and entropy, see Methods) was cal- 
culated (Borst and Theunissen, 1999; Sadeghi et al., 2007; Ince 
et al., 2009). The coding fraction of number of spikes, time win- 
dow, and time to first spike was relatively small (see Table 4) 
demonstrating a remarkable redundancy in the parameter space. 
Once inhibition was removed by blocking al and a 6 recep- 
tors, the coding fraction increased for time window and time 
to first spike, but not apparently so for the number of spikes 
(Table 4), indicating a preferential control of inhibition over spike 
patterning and timing rather than on spike number per se. 



INFORMATION CODING IN THE EXCITATORY/INHIBITORY SPACE 

In order to assess the impact of inhibition on information cod- 
ing (i.e., to what extent the same responses could be obtained 
with different parameters), the coding fraction (actually the ratio 



THE IMPACT OF E/l PHASE ON SPIKE GENERATION AND ON NMDA 
AND VDCC CURRENTS 

Beside the coordinated action of excitation and inhibition caused 
by concurrent mossy fiber activation of granule cells and Golgi 
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FIGURE 5 | The role of E/l and pE/pl in controlling granule cell firing 
parameters in response to mossy fiber bursts. (A) (Top) The relative 
changes of number of spikes (ns), time window (tw), time to first spike 
{ttfs) and first spike standard deviation [fssd) in control (con, black) and after 
switch off of a1 (a1-OFF, gray) or a6 (a6-OFF, light gray) receptors is shown 
for the low-£// and high- E/l quadrants identified in Figure 4B. It should be 
noted that the E/l balance primarily affects ns and tw.(Bottom) The relative 
change of number of spikes (ns), time window (tw), time to first spike (ttfs) 
and first spike standard deviation (fssd) is compared with phasic inhibition 
off (phasic OFF) and tonic inhibition off (tonic OFF). Blocking tonic inhibition 
does not influence notably the effect of phasic inhibition on the time-related 
spiking parameters, time to first spike and first spike standard deviation. 
(B) The relative changes of number of spikes (ns), time window (tw), time 
to first spike (ttfs) and first spike standard deviation (fssd) are reported as a 
function of the inhibitory (and excitatory release probability (Pi, Pe), 
showing that Pi affects modestly the parameters, while raising Pe 
remarkably increases number of spikes and decreases time to first spike 
and first spike standard deviation. 



cells shown in Figure 4, the excitatory and inhibitory synapses 
impinging on granule cells can work in any arbitrary phase rela- 
tionship. Given their different time-course and temporal integra- 
tion during trains (cf. Figures 2, 3), al and ot6 receptor-mediated 
currents might have differential impact on firing depending on 



the E/I phase (Figure 7A). Therefore, we have performed sim- 
ulations using different phase-lags of inhibition with respect 
to excitation generating tuning curves for the spiking parame- 
ters (in these simulations pE and pi were not changed since 
we used single stimulation pulses preventing their impact on 
short-term train dynamics). The most apparent changes occurred 
when a6 receptors were switched off. Interestingly, number of 
spikes, time window, time to first spike and first spike stan- 
dard deviation showed resonance around —10 ms phase-lag. 
When al receptors were switched off, resonance appeared around 
10 ms phase-lag. The impact of the E/I phase was evaluated 
also for currents generated by NMDA channels and voltage- 
dependent calcium channels (VDCCs) (Figure 7B). Again, the 
most apparent changes occurred when the a6 receptor was 
turned off. The VDCC current decreased monotonically with 
positive phase-lag while the NMDA current showed broad tun- 
ing distributed between —20 and +30 ms phase-lag. When the 
al receptor was turned off, resonance peaked around +10 ms 
phase-lag. 

These simulations showed that al and a6 receptor-mediated 
effects combined in tuning spike generation toward specific 
phase-lags and that a6 receptor-mediated effects were larger and 
broader than al receptor-mediated effects. 

MODULATION OF TUNING CURVES BY THE E/l PHASE AND TONIC 
INHIBITION 

The tuning curves reported in Figure 7 revealed dependence of 
granule cell responses on the E/I phase-lag. This effect was eval- 
uated by fitting the tuning curves with Lorentzian functions 
yielding Q = A/co, a parameter indicating the tendency of the 
response to peak at a certain phase-lag 9 max with area A and 
broadness co (Equation 4, see Methods; Figure 8A). 

Lorentzian parameters, which report the effectiveness of tun- 
ing toward a preferential phase-lag, showed specific properties 
(Figure 8B). For all parameters, Q(a6) > Q(al) reflected the 
larger impact of a6 than al receptor-mediated charge trans- 
fer, w(a6) > <w(al) reflected the longer duration of a6 than al 
receptor-mediated responses, and 0 max (a6) < 0 and f? mflx (al) > 
0 reflected the slower raise of a6 than al receptor-mediated 
responses (cf. Figures 2, 3). Interestingly, Q(a6) was much larger 
for time window, time to first spike and first spike standard 
deviation than for number of spikes (Figure 8B). This pattern was 
opposite to that observed for burst transmission (cf. Figure 5A), 
implying that only the time-dependent parameters were tuned by 
the E/I phase difference. It should also be noted that Q(a6) of the 
NMDA current was very small, reflecting the slow time course of 
this current and its broad integration times (DAngelo et al., 1995; 
Schwartz et al, 2012). 

Modulation of Q(a6) for time window, time to first spike and 
first spike standard deviation tended to decrease with E/I, accord- 
ing to the powerful effect of phasic inhibition in controlling time- 
dependent parameters (Figure 8C). Modulation of Q(a6) with 
E/I for the VDCC current followed strictly that of the number of 
spikes (DAngelo et al, 1999; Gandolfi et al., 2014). Modulation of 
Q(a6) with E/I for the NMDA current was marginal but, specifi- 
cally, 9-max (a6) shifted progressively from —20 to +30 ms (phase- 
lag precession, cf. Figure 7B) with increasing E/I (Figure 8D) 
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FIGURE 6 | Model predictions of experimental GABA-A receptor 
blockage. (A) The experimental parameters reported in Figures 1C,E are 

reported on the control E/I response space for the 5 cells showing E < 4 and 
/ < 4 yielding the corresponding Pe and Pi values. The position of symbols 
and their size correspond to the mean ± MSE of data values extrapolated 
from number of spikes (ns), time window (tw), time to first spike {ttfs) and 
first spike standard deviation (fssd) (same gray scale as in B). (B) The Pi and 
Pe values obtained in A and averaged across the different parameters (sn, tw, 



ttfs, and fssd) are reported as mean ± MSE for each of the 5 cells considered 
in (A). The release probabilities felt in a narrow range around the average 
values of Pe = 0.51 ± 0.06 and Pi = 0.42 ± 0.04 (dashed lines). (C) The plot 
compares the relative change of number of spikes (ns), time window (tw), 
time to first spike (ttfs) and first spike standard deviation (fssd) after gabazine 
perfusion in the 5 recordings shown in (A,B) (mean ± MSE) with that caused 
by switching off phasic and tonic inhibition in the model. The data points 
were fitted with a linear function yielding a slope of 0.95 (R 2 = 0.914, n = 4). 



Table 4 | Coding fraction of number of spikes [sn), time window ( tw) 
and time to first spike (ttfs) was relatively small demonstrating a 
remarkable redundancy in the parameter space. 





sn 


tw 


ttfs 


CTRL 


0.123 


0.144 


0.192 


a1 OFF 


0.128 


0.165 


0.211 


0(6OFF 


0.121 


0.181 


0.233 



Once inhibition was removed by blocking a 7 and a 6 subunit-containing recep- 
tors, the coding fraction increased for time window and time to first spike, but 
not apparently so for the number of spikes, indicating a preferential control of 
inhibition over spike patterning and timing rather than on spike number per se. 



unveiling a complex influence of a6 receptor-mediated inhibition 
on voltage-dependent NMDA channel unblock. 

Finally, in order to determine the effect of tonic inhibition on 
granule cell responses to dynamic patterns, tuning curves were 
generated with tonic inhibition blocked. The plot in Figure 8C 
shows the difference between the cases in which tonic inhibition 
was on or off. The tonic current biased the effect of inhibition pro- 
portionately to the E/I balance without altering the dynamics of 
inhibition substantially. Therefore, the model predicts that tonic 
inhibition simply scales the effect of phasic inhibition on firing 
parameters. 



DISCUSSION 

In this paper we have developed a detailed glomerular microcir- 
cuit model in order to assess how phasic and tonic components of 
GABA-A receptor-mediated inhibitory neurotransmission con- 
tribute to dynamic regulation of firing in cerebellar granule cells. 
The main prediction of this model is that phasic inhibition can 
dynamically regulate granule cell firing in response to temporally 
organized patterns of mossy fiber and Golgi cell activity. Phasic 
inhibition preferentially controlled the number of spikes during 
transmission of short bursts, but it showed improved control over 
time-related parameters (time to first spike, first spike standard 
deviation and time window) when the E/I phase was changed. In 
all conditions, the overall impact of a6 was larger than that of al 
subunit-containing receptors. However, al receptors controlled 
granule cell responses in a narrow ±10 ms band while 016 recep- 
tors showed broader ±50 ms tuning. Eventually, phasic inhibition 
fine-tuned the initiation and termination of granule cell firing 
during burst-burst transmission, as well as the probability of fir- 
ing at different E/I phases, on the millisecond-scale. The model 
also predicts that tonic inhibition can bias these effects without 
changing their nature substantially. These results imply that the 
specific organization of phasic inhibitory mechanisms can effec- 
tively regulate the dynamics of spike generation at the mossy 
fiber - granule cell relay of cerebellum without the intervention 
of tonic inhibition. 
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FIGURE 7 | Effect of phasic inhibition on granule cell response with 
variable E/l phase. (A) Exemplar simulations of inhibition occurring at 
different phases compared to excitation (E/l phase lag of —10, 5, 10 ms, as 
indicated in the insets) in control and after ot1 or ct6 receptor switch-off. (B) 
The plots show the changes of number of spikes (ns), time window (tw), 
time to first spike (ttfs), first spike standard deviation (fssd), VDCC currents 



and NMDA currents occurring at different E/l phase lags. The data represent 
the differences (a1 0 FF - control) and (a6 0 FF"Control) for simulations carried out 
with 3-4 inhibitory fibers and 1-4 excitatory fibers (gray coded as indicated in 
the inset). Each line is the interpolation to average data obtained with 
different release probabilities. It should be noted that the data distributions 
show resonance around 0 ms. 



THE PREDICTIVE POWER OF THE GLOMERULAR MICR0CIRCUIT MODEL 

The glomerular microcircuit model could fit and explain the fir- 
ing patterns observed in experimental data following mossy fiber 
burst stimulation. Once experimental granule cell firing parame- 
ter (time to first spike, first spike standard deviation, number of 
spikes and time window) values were reported on the simulated 
response space, the glomerular microcircuit model yielded pE and 
pi values corresponding to those observed experimentally (aver- 
age pE = 0.51 vs. 0.55 in (Sola et al, 2004); average pi = 0.42 vs. 
0.42 in Mapelli et al., 2009). Moreover, the model properly fit- 
ted the spiking parameters changes caused by gabazine, a drug 
that blocks both phasic and tonic inhibition. Therefore the model 
was endowed with mechanisms capable of capturing the impact 
of inhibitory neurotransmission on granule cell firing parameters 
and could be used to predict the effect of different E/I phase lags. 

It should be noted that in our model the effect of metabotropic 
receptors was excluded to isolate the ionotropic component 



of inhibition. Experimentally, the application of mGluRs and 
GABA-B receptor blockers caused a net dis-inhibitory effect (cf. 
Figure 1), which could reflect multiple mechanisms. The mGluRs 
on presynaptic Golgi cell terminals might have reduced GABA 
release, and this would suffice to explain the increased number 
of spikes and duration of discharge in most of the recordings. 
Conversely, GABA-B receptor blockage on granule cells and on 
Golgi cell terminals (Rossi et al, 2006; Mapelli et al., 2009) would 
reduce granule cell excitability and enhance inhibitory neuro- 
transmission, thereby sorting the opposite effect. Future models 
of glomerula neurotransmission may be designed to include 
metabotropic receptor control (Mapelli et al., 2014). 

PROPERTIES OF THE INHIBITORY NEUROTRANSMISSION MODEL 

The inhibitory neurotransmission model was based on GABA- 
A receptor kinetics and was fitted to voltage-clamp recordings. 
Simulations of glomerular inhibition were consistent with a 
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FIGURE 8 | Resonance of the parameter distributions at different E/l 
levels. (A) Exemplar fitting with a Lorentzian function (Equation 4) applied to 
the phase-lag plot (cf. Figure 7B) indicating the geometrical meaning of the 
parameters Q, A, w, 9 max . (B) The average values of Lorentzian parameters 
during ot1 or a6 switch-off are reported for number of spikes (ns), time window 
(tw), time to first spike (ttfs), first spike standard deviation (fssd), VDCC current 
and NMDA current. Note that Q(u6) and <u(a6) is always larger than Q(a1 ) and 
n)(a1 ), while 8 max (a6) is negative and t9 msx (a1 ) is positive. These data are 



averaged over the whole parameter response space. (C) The top panel shows 
the Q factors of the aQ 0FF difference phase-lag curves (a6 0 FF"CTRL). Note that 
Qof time window (tw), time to first spike (ttfs) and first spike standard deviation 
(fssd) decrease with increasing E/l, while the other parameters are much less 
affected. The bottom panel shows the difference Q factors obtained as Q (tonic 
OFF)-Q (tonic ON). Note that the trends of these Q plots is qualitatively similar 
to that reported above. (D) Precession of 0 max (aS) for the NMDA Q factor. The 
phase of the peak (color coded) increases with increasing E/l. 



mechanism in which oti receptors determine the rapid response 
to released GABA characterizing the direct component, while ag 
receptors determine the slow temporal integration of spillover 
currents characterizing the indirect component (Cherubini and 
Conti, 2001; Farrant and Nusser, 2005; Glykys and Mody, 2007). 
The model predicts that, while oil receptor-mediated responses 
can precisely follow GABA concentration changes, c*6 receptor- 
mediated responses tend to saturate early during the train. Thus, 
as observed experimentally, raise in p determined minor changes 
in synaptic charge transfer matching experimental observations 
(Mapelli et al, 2009). 

The adherence of the inhibitory neurotransmission model 
to the microscopic properties of the synapse could be defined 



at different levels. First, the GABA-A receptor kinetic schemes 
(derived from Pugh and Raman, 2005) appropriately reproduced 
elPSCs and their changes in different functional conditions, sug- 
gesting that GABA affinities and the balance between open, closed 
and desensitized states were in the physiological range for granule 
cell receptors. Secondly, by setting the same diffusion coefficient 
for GABA and glutamate (DiGregorio et al., 2002), the number 
of diffusing molecules and the diffusion distance were similar to 
those estimated for the mossy fiber — granule cell synapse (Nieus 
et al., 2006). This was not unexpected, given the substantial sim- 
ilarity in the number of releasing sites involved and the common 
diffusion space for the two neurotransmitters. Thirdly, the p val- 
ues obtained by fitting the model to experimental data were 
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almost identical to those estimated experimentally with quan- 
tal analysis and no changes in the rate of vesicle recycling or in 
any other model parameters were needed to fit the experimental 
elPSC trains. 

INHIBITORY CONTROL OF GRANULE CELL FIRING DURING BURST 
TRANSMISSION 

Model simulations explained the complex effect of phasic inhi- 
bition on the burst-burst transmission modality of the granular 
layer (Chadderton et al, 2004; Rancz et al, 2007; Arenz et al., 
2008). The E/I balance was most effective in controlling granule 
cell tonic discharge (the ~9-fold number of spikes increase was 
associated with a ~3-fold time window increase, so that the out- 
put spike frequency increased by about 3 times) but had minor 
effects on time-related parameters, time to first spike and first 
spike standard deviation. The pE/pI balance was most effective 
in controlling number of spikes, time to first spike and first spike 
standard deviation. Raising pE specifically increased number of 
spikes and reduced time to first spike, consistent with the fact 
that AMPA and NMDA receptor activation shows a remarkable 
modulation with pi: (Nieus et al., 2006). Raising pi depressed all 
parameters but its effect was small, probably because inhibition 
during a burst largely depends on the charge carried by ot6 recep- 
tors, which are early saturated (see above and cf. Figures 2, 3). 
The first spike standard deviation was almost equally reduced by 
raising pE and pi, probably reflecting reduction of synaptic noise 
determined by the larger number of released quanta (Sola et al., 
2004; Mapelli et al, 2009). It turns out that, during burst trans- 
mission, the major effect of inhibition is to control the number 
and frequency of spikes emitted by granule cells through a reg- 
ulation of the relative number of active Golgi cell synapses. Both 
E/I and pE/pI effects were enhanced by GABA receptor switch-off, 
but the contribution of ot6 was always much larger than that of al 
receptors reflecting the fact that a6 receptors controlled most of 
the inhibitory charge transfer. 

The tonic GABA leakage in the model did not contribute to 
the effect of inhibitory transmission substantially beside propor- 
tionately increasing number of spikes and time window. The tonic 
GABA leakage was modeled on the basis of experimental obser- 
vations in cerebellar slices. The cerebellar glomerulus is an ideal 
structure to preserve tonic inhibition in slices provided that the 
glial sheath is intact (Hamann et al., 2002; Mitchell and Silver, 
2003; Mapelli et al, 2014). Although the topic is still debated, 
tonic inhibition was reported in vivo revealing properties similar 
as in slices (Chadderton et al., 2004), suggesting that predic- 
tions made by our model could apply to synaptic transmission 
in natural conditions. 

INHIBITORY CONTROL OF GRANULE CELL FIRING WITH VARIABLE E/I 
PHASE 

Model simulations predicted a specific role for phasic inhibition 
in controlling granule cell responses with different E/I phases 
lags. The most apparent regulatory effects on granule cell firing 
occurred when the u.6 receptor-mediated current was turned off. 
This was due to the long time course and large charge trans- 
fer of a6 receptor-mediated responses, extending the influence 
of inhibition for about ±50 ms around granule cell excitation. 



When the al receptor-mediated current was turned off, the effects 
were smaller and resonance was restricted over about ±10 ms. 
Interestingly, time-related firing parameters (time window, time 
to first spike, first spike standard deviation) showed much sharper 
tuning than tonic firing (number of spikes). Moreover, tuning 
of time-related parameters decreased with E/I along with the 
decreasing inhibitory control, while GABA leakage did not alter 
the dynamics of inhibition substantially. 

The model also predicted how VDCC and NMDA currents 
would be regulated by al and a6 receptors. VDCC current regula- 
tion strictly followed the number of spikes, as these channels open 
during the upstroke (D'Angelo et al., 1997, 2001; Gandolfi et al., 
2014). The a6 receptor-mediated regulation of the NMDA cur- 
rent showed a broad tuning extending over ±50 ms and reflecting 
the slow NMDA current kinetics (D'Angelo et al., 1995; Schwartz 
et al, 2012). Puzzlingly, the effect of inhibition on the NMDA 
current showed a precession with E/I (from —20 to ±30 ms 
phase-lag) suggesting a complex voltage-dependent interaction 
between NMDA channel unblock, depolarization and inhibitory 
control. 

IMPLICATIONS FOR DYNAMIC INHIBITORY CONTROL DURING CIRCUIT 
ACTIVITY 

The dynamic control of inhibition by a6 and al receptors could 
have specific consequences during local microcircuit activation. 
Feed-forward inhibition affects granule cell excitation around 
5 ms after mossy fiber activation (D'Angelo and De Zeeuw, 2009). 
With 5 ms phase-lag, the model predicted a higher impact for a6 
mediated inhibition (except for the number of spikes at low E/I, 
where a6 and al affected the response in a similar way). Feed- 
back inhibition occurs after about 10 ms. With 10 ms phase-lag, 
the overall effect was also dominated by a6 mediated inhibi- 
tion, but al and a6 effects were similar for time-to-first spike 
and first-spike standard deviation (and also on the number of 
spikes at low E/I). Therefore, in summary, the a6 component 
dominated on feed-forward inhibition, while the al component 
gained relevance in controlling the timing of granule cell output 
firing at low E/I during feed-back inhibition. Interestingly, during 
feedback inhibition at low E/I the al component is often more 
effective than the u6 component. Concerning the NMDA current, 
phase-lag precession suggested that the maximum impact of a6 
mediated inhibition at 5-10 ms phase-lag was also achieved at low 
E/I. The multiple combinations of these effects and the complex 
inhibitory patterns provided by al and a6 subunit-containing 
GABA-A receptors make phasic inhibition particularly suitable to 
finely regulate the delay and precision of granule cell firing. 

These properties may have an impact on spike-timing depen- 
dent plasticity (STDP) induction and on its interaction with 
low-frequency coherent oscillations reported in the cerebellar 
granular layer (Pellerin and Lamarre, 1997; Hartmann and Bower, 
1998) and deserve future experimental assessment. 

CONSIDERATIONS ON INHIBITORY CONTROL OF INFORMATION 
TRANSFER 

The extensive simulation strategy adopted in this paper suggests 
how phasic inhibition provided by Golgi cells controls informa- 
tion transfer through the granular layer. For each E/I value, release 
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probabilities generated a complex subspace, such that similar fir- 
ing parameter values could be found with different E/I and pE/pI 
combinations with a redundancy of 80 -r- 90% (CF = 0.1 -f- 0.2) 
(Borst and Theunissen, 1999; Sadeghi et al., 2007; Ince et al., 
2009). The actual number of excitatory synapses may be relatively 
unimportant, provided that the number of inhibitory synapses 
is properly set, suggesting that the granular layer has intrinsic 
re-scaling properties with respect to the intensity of the inputs. 
Synaptic plasticity, mostly through pE, could allow to select points 
in this computational landscape for any specific E/I and pE/pI 
combination. 

The amount of transmitted information (MI; Borst and 
Theunissen, 1999; Sadeghi et al, 2007; Ince et al, 2009), depends 
on the signal-to-noise (S/N) ratio. Since short bursts (1-3 spikes: 
see Figure 4B) at the mossy fiber — granule cell synapse carry 
the largest fraction of MI (Arleo et al., 2010), it is expected 
that regulation of first-spike delay and precision by phasic inhi- 
bition contributes substantially to regulate MI transfer through 
the cerebellum input stage (Garrido et al., 2013). It should also 
be noted that quantal fluctuations in GABA release generate 
synaptic noise, which should be smaller the larger the num- 
ber of active synapses and the higher their release probabilities. 
Therefore, redundant E/I and pE/pI combinations giving same 
granule cell output patterns may not be equivalent in terms of MI 
transfer. 

THE FIVE REASONS WHY PHASIC INHIBITION IS NEEDED FOR 
DYNAMIC CONTROL 

In summary, simulations identify five main reasons why phasic 
inhibition is needed for dynamic control of spike transmission at 
the mossy fiber — granule cell relay. (1) During burst transmis- 
sion, tonic inhibition simply scales the effect of phasic inhibition 
on tonic spike discharge (number of spikes) without affecting 
time-related parameters (time to first spike and time window) 
significantly. (2) During inhibition with variable E/I phase, the 
pattern observed during bursts is reversed, so that phasic inhibi- 
tion specifically affects time-related parameters without affecting 
tonic spike discharge significantly. (3) During inhibition with 
variable E/I phase, tonic GABA leakage does not alter the dynam- 
ics of inhibition substantially, so that tonic inhibition simply 
scales the effect of phasic inhibition on firing parameters. (4) 
The differential kinetics imposed by al and a6 receptors dur- 
ing inhibition with bursts or with variable E/I phase may not 
be explained by tonic inhibitory effects. (5) Phasic inhibition 
reduced redundancy specifically for time-related parameters. 

CONCLUSIONS 

Biologically realistic simulations show that phasic inhibition, by 
exploiting the different kinetics of al and a6 receptor-mediated 
responses, can effectively regulate the delay and precision of 
granule cell firing. The effect of a6 receptors dominated during 
burst transmission while that of al receptors emerged during 
single-pulse transmission with variable phase lag, so that al and 
a6 receptors proved suitable for single-spike and burst control, 
respectively. These mechanisms in turn may regulate coincidence 
detection of calcium signals required for STDP and for its inte- 
gration during theta-cycles (DAngelo and De Zeeuw, 2009). 



Conversely, no specific effects of tonic inhibition occurred beyond 
a proportionate scaling of granule cell discharge duration and fre- 
quency, compatible with a gain control mechanisms over tonic 
firing (Mitchell and Silver, 2003). Therefore, these model simu- 
lations imply that phasic inhibition endows the cerebellar cortex 
with a mechanism suitable for improving spatio-temporal recon- 
figuration of the mossy fiber input and for promoting dynamic 
signal processing in the cerebellar cortex (Marr, 1969; Eccles, 
1973; Fujita, 1982). Critical tests to this prediction could be 
obtained by experiments, in which multiple different mossy fiber 
and Golgi cell axon fibers are activated independently with arbi- 
trary phases. Although this protocol would be impractical with 
current electrophysiology, it may become at hand by combining 
recently developed imaging and optogenetic techniques (e.g., see 
Gandolfi et al., 2014). 

ACKNOWLEDGMENTS 

We thank Leda Roggeri for technical assistance. This work 
was supported by European Union grants to Egidio DAngelo 
[CEREBNET FP7-ITN238686, REALNET FP7-ICT270434, 
Human Brain Project (HBP-604102)]. 

REFERENCES 

Arenz, A., Silver, R. A., Schaefer, A. T., and Margrie, T. W. (2008). The contribution 
of single synapses to sensory representation in vivo. Science 321, 977-980. doi: 
10.1126/science.ll58391 

Arleo, A., Nieus, T., Bezzi, M., D'Errico, A., D'Angelo, E., and Coenen, O. J. (2010). 
How synaptic release probability shapes neuronal transmission: information- 
theoretic analysis in a cerebellar granule cell. Neural Comput. 22, 2031-2058. 
doi: 10.1162/NECO_a_00006-Arleo 

Armano, S., Rossi, R, Taglietti, V., and D'Angelo, E. (2000). Long-term potentiation 
of intrinsic excitability at the mossy fiber-granule cell synapse of rat cerebellum. 
/. Neurosci. 20, 5208-5216. 

Barberis, A., Cherubini, E., and Mozrzymas, J. W. (2000). Zinc inhibits minia- 
ture GABAergic currents by allosteric modulation of GABAA receptor gating. 
/. Neurosci. 20, 8618-8627. 

Borst, A., and Theunissen, F. E. (1999). Information theory and neural coding. Nat. 
Neurosci. 2, 947-957. doi: 10.1038/14731 

Brandalise, E, Mapelli, L., Gerber, U., and Rossi, P. (2012). Golgi cell-mediated 
activation of postsynaptic GABA(B) receptors induces disinhibition of the 
Golgi cell-granule cell synapse in rat cerebellum. PLoS ONE 7:e43417. doi: 
10.1371/journal.pone.0043417 

Brickley, S. G., Cull-Candy, S. G., and Farrant, M. (1996). Development of 
a tonic form of synaptic inhibition in rat cerebellar granule cells result- 
ing from persistent activation of GABAA receptors. /. Physiol. 497(Pt 3), 
753-759. 

Brickley, S. G., ReviJla, V., Cull-Candy, S. G., Wisden, W., and Farrant, M. (2001). 
Adaptive regulation of neuronal excitability by a voltage-independent potas- 
sium conductance. Nature 409, 88-92. doi: 10.1038/35051086 

Cesana, E., Pietrajtis, K., Bidoret, C, Isope, P., D'Angelo, E., Dieudonne, S., et al. 
(2013). Granule cell ascending axon excitatory synapses onto Golgi cells imple- 
ment a potent feedback circuit in the cerebellar granular layer. /. Neurosci. 33, 
12430-12446. doi: 10.1523/JNEUROSCI.4897-11.2013 

Chadderton, P., Margrie, T. W., and Hausser, M. (2004). Integration of quanta in 
cerebellar granule cells during sensory processing. Nature 428, 856-860. doi: 
10.1038/nature02442 

Cherubini, E., and Conti, F. (2001). Generating diversity at GABAergic synapses. 
Trends Neurosci. 24, 155-162. doi: 10.1016/S0166-2236(00)01724-0 

D'Angelo, E., De Filippi, G., Rossi, P., and Taglietti, V. (1995). Synaptic excitation 
of individual rat cerebellar granule cells in situ: evidence for the role of NMDA 
receptors. /. Physiol. 484(Pt 2), 397^13. 

D'Angelo, E., De Filippi, G., Rossi, P., and Taglietti, V. (1997). Synaptic activa- 
tion of Ca2-|- action potentials in immature rat cerebellar granule cells in situ. 
J. Neurophysiol. 78, 1631-1642. 



Frontiers in Cellular Neuroscience 



www.frontiersin.org 



August 2014 | Volume 8 | Article 246 | 16 



Nieus et al. 



Phasic inhibition in cerebellar granule cells 



D'Angelo, E., and De Zeeuw, C. I. (2009). Timing and plasticity in the cerebellum: 
focus on the granular layer. Trends Neurosci. 32, 30-40. doi: 10.1016/j.tins.2008. 
09.007 

D'Angelo, E„ Nieus, T., Maffei, A., Armano, S., Rossi, P., Taglietti, V., et al. (2001). 
Theta-frequency bursting and resonance in cerebellar granule cells: experimen- 
tal evidence and modeling of a slow k + -dependent mechanism. /. Neurosci. 21, 
759-770. 

DAngelo, E., Rossi, P., Armano, S., and Taglietti, V. (1999). Evidence for NMDA 
and mGlu receptor-dependent long-term potentiation of mossy fiber-granule 
cell transmission in rat cerebellum. /. Neurophysiol. 81, 277-287. 

DAngelo, E., Rossi, P., and Taglietti, V. (1993). Different proportions of N- 
methyl-D-aspartate and non-iV-methyl-D-aspartate receptor currents at the 
mossy fibre-granule cell synapse of developing rat cerebellum. Neuroscience 53, 
121-130. doi: 10.1016/0306-4522(93)90290-V 

DAngelo, E., Solinas, S., Mapelli, J., Gandolfi, D., Mapelli, L., and Prestori, F. 

(2013) . The cerebellar Golgi cell and spatiotemporal organization of granular 
layer activity. Front. Neural Circuits 7:93. doi: 10.3389/fncir.2013.00093 

D'Errico, A., Prestori, F., and DAngelo, E. (2009). Differential induction of bidi- 
rectional long-term changes in neurotransmitter release by frequency-coded 
patterns at the cerebellar input. /. Physiol. 587, 5843-5857. doi: 10.1113/jphys- 
iol.2009.177162 

Dieudonne, S. (1998). Submillisecond kinetics and low efficacy of parallel fibre- 

Golgi cell synaptic currents in the rat cerebellum. /. Physiol. 510(Pt 3), 845-866. 

doi: 10.1 1 1 1/j. 1469-7793. 1998.845bj.x 
DiGregorio, D. A., Nusser, Z., and Silver, R. A. (2002). Spillover of glutamate onto 

synaptic AMPA receptors enhances fast transmission at a cerebellar synapse. 

Neuron 35, 521-533. doi: 10.1016/S0896-6273(02)00787-0 
Diwakar, S., Magistretti, J., Goldfarb, M., Naldi, G, and DAngelo, E. (2009). 

Axonal Na + channels ensure fast spike activation and back-propagation in 

cerebellar granule cells. /. Neurophysiol. 101, 519-532. doi: 10.1152/jn.90382. 

2008 

Dugue, G. P., Brunei, N., Hakim, V., Schwartz, E., Chat, M., Levesque, M., et al. 
(2009). Electrical coupling mediates tunable low-frequency oscillations and res- 
onance in the cerebellar Golgi cell network. Neuron 61, 126-139. doi: 10.1016/j. 
neuron.2008.1 1.028 

Eccles, J. C. (1973). The cerebellum as a computer: patterns in space and time. 
/. Physiol. 229, 1-32. 

Eccles, J., Llinas, R., and Sasaki, K. (1964). Golgi cell inhibition in the cerebellar 

cortex. Nature 204, 1265-1266. doi: 10.1038/2041265a0 
Farrant, M., and Nusser, Z. (2005). Variations on an inhibitory theme: phasic and 

tonic activation of GABA(A) receptors. Nat. Rev. Neurosci. 6, 215-229. doi: 

10.1038/nrnl625 

Forti, L., Cesana, E., Mapelli, J., and D'Angelo, E. (2006). Ionic mechanisms of 
autorhythmic firing in rat cerebellar Golgi cells. /. Physiol. 574, 711-729. doi: 
10.1 1 13/jphysiol.2006. 1 10858 

Fujita, M. (1982). Adaptive filter model of the cerebellum. Biol. Cybern. 45, 
195-206. doi: 10.1007/BF00336192 

Gandolfi, D., Lombardo, P., Mapelli, J., Solinas, S., and D'Angelo, E. (2013). Theta- 
frequency resonance at the cerebellum input stage improves spike timing on 
the millisecond time-scale. Front. Neural Circuits 7:64. doi: 10.3389/fncir.2013. 
00064 

Gandolfi, D., Pozzi, P., Tognolina, M., Chirico, G., Mapelli, J., and D'Angelo, E. 

(2014) . The spatiotemporal organization of cerebellar network activity resolved 
by two-photon imaging of multiple single neurons. Front. Cell. Neurosci. 8:92. 
doi: 10.3389/fncel.2014.00092 

Garrido, J. A., Ros, E., and D'Angelo, E. (2013). Spike timing regulation on the mil- 
lisecond scale by distributed synaptic plasticity at the cerebellum input stage: a 
simulation study. Front. Comput. Neurosci. 7:64. doi: 10.3389/fncom.2013.00064 

Glykys, J., and Mody, I. (2007). Activation of GABAA receptors: views from 
outside the synaptic cleft. Neuron 56, 763-770. doi: 10.1016/j.neuron.2007. 
11.002 

Griffen, T. C., and Maffei, A. (2014). GABAergic synapses: their plasticity and role 

in sensory cortex. Front. Cell. Neurosci. 8:91. doi: 10.3389/fncel.2014.00091 
Hadley, S. H., and Amin, J. (2007). Rat alpha6beta2delta GABAA receptors exhibit 

two distinct and separable agonist affinities. /. Physiol. 581, 1001-1018. doi: 

10.1 1 13/jphysiol.2007.132886 
Hamann, M., Rossi, D. J., and Attwell, D. (2002). Tonic and spillover inhibition 

of granule cells control information flow through cerebellar cortex. Neuron 33, 

625-633. doi: 10.1016/S0896-6273(02)00593-7 



Hamori, J., and Somogyi, J. (1983). Differentiation of cerebellar mossy fiber 
synapses in the rat: a quantitative electron microscope study. /. Comp. Neurol. 
220, 365-377. doi: 10.1002/cne.902200402 

Hartmann, M. J., and Bower, J. M. (1998). Oscillatory activity in the cerebellar 
hemispheres of unrestrained rats. /. Neurophysiol. 80, 1598-1604. 

Hines, M. L., and Carnevale, N. T. (2001). NEURON: a tool for neuroscientists. 
Neuroscientist 7, 123-135. doi: 10.1177/107385840100700207 

Hull, C, and Regehr, W. G. (2012). Identification of an inhibitory circuit that reg- 
ulates cerebellar Golgi cell activity. Neuron 73, 149-158. doi: 10.1016/j.neuron. 
2011.10.030 

Ince, R. A., Petersen, R. S., Swan, D. C, and Panzeri, S. (2009). Python for 

information theoretic analysis of neural data. Front. Neuroinform. 3:4. doi: 

10.3389/neuro.ll.004.2009 
Jakab, R. L., and Hamori, J. (1988). Quantitative morphology and synaptology of 

cerebellar glomeruli in the rat. Anat. Embryol. (Bed.) 179, 81-88. doi: 10.1007/ 

BF00305102 

Jones, M. V., Jonas, P., Sahara, Y., and Westbrook, G. L. (2001). Microscopic kinet- 
ics and energetics distinguish GABA(A) receptor agonists from antagonists. 
Biophys. J. 81, 2660-2670. doi: 10.1016/S0006-3495(01)75909-7 

Kanichay, R. T, and Silver, R. A. (2008). Synaptic and cellular properties of the 
feedforward inhibitory circuit within the input layer of the cerebellar cortex. 
/. Neurosci. 28, 8955-8967. doi: 10.1523/JNEUROSCI.5469-07.2008 

Lamsa, K. P., Kullmann, D. M., and Woodin, M. A. (2010). Spike-timing 
dependent plasticity in inhibitory circuits. Front. Synaptic Neurosci. 2:8. doi: 
10.3389/fnsyn.2010.00008 

Llinas, R., Urbano, F. )., Leznik, E., Ramirez, R. R., and van Marie, H. J. (2005). 
Rhythmic and dysrhythmic thalamocortical dynamics: GABA systems and the 
edge effect. Trends Neurosci. 28, 325-333. doi: 10.1016/j.tins.2005.04.006 

Maex, R., and De Schutter, E. (1998). Synchronization of golgi and granule 
cell firing in a detailed network model of the cerebellar granule cell layer. 
/. Neurophysiol. 80, 2521-2537. 

Mann, E. O., and Paulsen, O. (2007). Role of GABAergic inhibition in hippocam- 
pal network oscillations. Trends Neurosci. 30, 343-349. doi: 10.1016/j. tins. 2007. 
05.003 

Mapelli, J., and D'Angelo, E. (2007). The spatial organization of long-term synap- 
tic plasticity at the input stage of cerebellum. /. Neurosci. 27, 1285-1296. doi: 
10.1 523/JNEUROSCI.4873-06.2007 

Mapelli, J., Gandolfi, D., and D'Angelo, E. (2010a). Combinatorial responses con- 
trolled by synaptic inhibition in the cerebellum granular layer. /. Neurophysiol. 
103, 250-261. doi: 10.1152/jn.00642.2009 

Mapelli, J., Gandolfi, D., and D'Angelo, E. (2010b). High-pass filtering and dynamic 
gain regulation enhance vertical bursts transmission along the mossy fiber 
pathway of cerebellum. Front. Cell Neurosci. 4:14. doi: 10.3389/fncel.2010. 
00014 

Mapelli, L., Rossi, P., Nieus, T., and D'Angelo, E. (2009). Tonic activation of GABAB 
receptors reduces release probability at inhibitory connections in the cerebellar 
glomerulus./. Neurophysiol. 101, 3089-3099. doi: 10.1152/jn.91190.2008 

Mapelli, L., Solinas, S., and D'Angelo, E. (2014). Integration and regulation 
of glomerular inhibition in the cerebellar granular layer circuit. Front. Cell. 
Neurosci. 8:55. doi: 10.3389/fncel.2014.00055 

Marr, D. (1969). A theory of cerebellar cortex. /. Physiol. 202, 437-470. 

Mitchell, S. J., and Silver, R. A. (2003). Shunting inhibition modulates neuronal 
gain during synaptic excitation. Neuron 38, 433-445. doi: 10.1016/S0896- 
6273(03)00200-9 

Nielsen, T. A., DiGregorio, D. A., and Silver, R. A. (2004). Modulation of glutamate 
mobility reveals the mechanism underlying slow-rising AMPAR EPSCs and the 
diffusion coefficient in the synaptic cleft. Neuron 42, 757-771. doi: 10.1016/j. 
neuron.2004.04.003 

Nieus, T, Sola, E., Mapelli, J., Saftenku, E., Rossi, P., and D'Angelo, E. (2006). LTP 
regulates burst initiation and frequency at mossy fiber-granule cell synapses 
of rat cerebellum: experimental observations and theoretical predictions. 
/. Neurophysiol. 95, 686-699. doi: 10.1152/jn.00696.2005 

Nusser, Z., and Mody, I. (2002). Selective modulation of tonic and phasic inhibi- 
tions in dentate gyrus granule cells. /. Neurophysiol. 87, 2624-2628. doi: 10.1152/ 
jn.00866.2001 

Nusser, Z., Roberts, J. D., Baude, A., Richards, J. G, and Somogyi, P. (1995). Relative 
densities of synaptic and extrasynaptic GABAA receptors on cerebellar gran- 
ule cells as determined by a quantitative immunogold method. /. Neurosci. 15, 
2948-2960. 



Frontiers in Cellular Neuroscience 



www.frontiersin.org 



August 2014 | Volume 8 | Article 246 | 17 



Nieus et al. 



Phasic inhibition in cerebellar granule cells 



Nusser, Z., Sieghart, W., and Somogyi, P. (1998). Segregation of different GABAA 

receptors to synaptic and extrasynaptic membranes of cerebellar granule cells. 

/. Neurosci. 18, 1693-1703. 
Pellerin, J. P., and Lamarre, Y. (1997). Local field potential oscillations in 

primate cerebellar cortex during voluntary movement. /. Neurophysiol. 78, 

3502-3507. 

Petrini, E. M., Nieus, T., Ravasenga, T., Succol, F., Guazzi, S., Benfenati, R, et al. 
(2011). Influence of GABAAR monoliganded states on GABAergic responses. 
/. Neurosci. 31, 1752-1761. doi: 10.1523/JNEUROSCI.1453-10.2011 

Pugh, J. R., and Raman, I. M. (2005). GABAA receptor kinetics in the cerebellar 
nuclei: evidence for detection of transmitter from distant release sites. Biophys. 
J. 88, 1740-1754. doi: 10.1529/biophysj. 104.055814 

Rancz, E. A., Ishikawa, T., Duguid, I., Chadderton, P., Mahon, S., and Hausser, 
M. (2007). High-fidelity transmission of sensory information by single cere- 
bellar mossy fibre boutons. Nature 450, 1245-1248. doi: 10.1038/nature 
05995 

Rossi, D. J., and Hamann, M. (1998). Spillover-mediated transmission at inhibitory 
synapses promoted by high affinity alpha6 subunit GABA(A) receptors 
and glomerular geometry. Neuron 20, 783-795. doi: 10.1016/S0896-6273(00) 
81016-8 

Rossi, D. J., Hamann, M., and Attwell, D. (2003). Multiple modes of GABAergic 
inhibition of rat cerebellar granule cells. /. Physiol. 548, 97-110. doi: 
10. 1 1 13/jphysiol.2002.036459 

Rossi, P., DAngelo, E., Magistretti, J., Toselli, M., and Taglietti, V. (1994). Age- 
dependent expression of high-voltage activated calcium currents during cere- 
bellar granule cell development in situ. Pflugers Arch. 429, 107-116. doi: 
10.1007/BF02584036 

Rossi, P., Mapelli, L., Roggeri, L., Gall, D., de Kerchove d'Exaerde, A., Schiffmann, S. 
N., et al. (2006). Inhibition of constitutive inward rectifier currents in cerebellar 
granule cells by pharmacological and synaptic activation of GABA receptors. 
Eur. J. Neurosci. 24, 419-432. doi: 10.1111/j.l460-9568.2006.04914.x 

Sadeghi, S. G, Chacron, M. J., Taylor, M. C, and Cullen, K. E. (2007). 
Neural variability, detection thresholds, and information transmission in the 
vestibular system. /. Neurosci. 27, 771-781. doi: 10.1523/JNEUROSCI.4690- 
06.2007 

Saviane, C, and Silver, R. A. (2007). Estimation of quantal parameters with 

multiple-probability fluctuation analysis. Methods Mol. Biol. 403, 303-317. doi: 

10.1007/978-l-59745-529-9_19 
Schwartz, E. J., Rothman, J. S., Dugue, G. P., Diana, M., Rousseau, C, Silver, 

R. A., et al. (2012). NMDA receptors with incomplete Mg 2+ block enable 

low-frequency transmission through the cerebellar cortex. /. Neurosci. 32, 

6878-6893. doi: 10.1523/JNEUROSCI.5736-11.2012 
Silver, R. A., Traynelis, S. E, and Cull-Candy, S. G. (1992). Rapid-time-course 

miniature and evoked excitatory currents at cerebellar synapses in situ. Nature 

355, 163-166. doi: 10.1038/355163a0 



Sola, E., Prestori, E, Rossi, P., Taglietti, V., and DAngelo, E. (2004). Increased 
neurotransmitter release during long-term potentiation at mossy fibre-granule 
cell synapses in rat cerebellum. /. Physiol. 557, 843-861. doi: 10.1113/jphys- 
iol.2003.060285 

Solinas, S., Forti, L., Cesana, E., Mapelli, J., De Schutter, E., and DAngelo, E. 
(2007a). Computational reconstruction of pacemaking and intrinsic electrore- 
sponsiveness in cerebellar Golgi cells. Front. Cell. Neurosci. 1:2. doi: 10.3389/ 
neuro.03.002.2007 

Solinas, S., Forti, L., Cesana, E., Mapelli, J., De Schutter, E., and DAngelo, E. 

(2007b). Fast-reset of pacemaking and theta-frequency resonance patterns in 

cerebellar golgi cells: simulations of their impact in vivo. Front. Cell. Neurosci. 

1:4. doi: 10.3389/neuro.03.004.2007 
Solinas, S., Nieus, T., and DAngelo, E. (2010). A realistic large-scale model of the 

cerebellum granular layer predicts circuit spatio-temporal filtering properties. 

Front. Cell. Neurosci. 4:12. doi: 10.3389/fncel.2010.00012 
Tia, S., Wang, J. E, Kotchabhakdi, N., and Vicini, S. (1996). Distinct deac- 
tivation and desensitization kinetics of recombinant GABAA receptors. 

Neuropharmacology 35, 1375-1382. doi: 10.1016/S0028-3908(96)00018-4 
Timmann, D., Watts, S., and Hore, J. (1999). Failure of cerebellar patients to time 

finger opening precisely causes ball high-low inaccuracy in overarm throws. 

/. Neurophysiol. 82, 103-114. 
Vos, B. P., Volny-Luraghi, A., and De Schutter, E. (1999). Cerebellar Golgi cells 

in the rat: receptive fields and timing of responses to facial stimulation. Eur. 

J. Neurosci. 11, 2621-2634. doi: 10.1046/j.l460-9568.1999.00678.x 
Wall, M. J. (2002). Furosemide reveals heterogeneous GABA(A) receptor expression 

at adult rat Golgi cell to granule cell synapses. Neuropharmacology 43, 737-749. 

doi: 10.101 6/S0028-3908(02)00085-0 
Wenner, P. (2011). Mechanisms of GABAergic homeostatic plasticity. Neural Plast. 

2011, 489470. doi: 10.1155/2011/489470 

Conflict of Interest Statement: The authors declare that the research was con- 
ducted in the absence of any commercial or financial relationships that could be 
construed as a potential conflict of interest. 

Received: 18 June 2014; accepted: 04 August 2014; published online: 25 August 2014. 
Citation: Nieus TR, Mapelli L and DAngelo E (2014) Regulation of output spike pat- 
terns by phasic inhibition in cerebellar granule cells. Front. Cell. Neurosci. 8:246. doi: 
10.3389/fncel.2014.00246 

This article was submitted to the journal Frontiers in Cellular Neuroscience. 
Copyright © 2014 Nieus, Mapelli and DAngelo. This is an open-access article dis- 
tributed under the terms of the Creative Commons Attribution License (CC BY). The 
use, distribution or reproduction in other forums is permitted, provided the original 
author(s) or licensor are credited and that the original publication in this jour- 
nal is cited, in accordance with accepted academic practice. No use, distribution or 
reproduction is permitted which does not comply with these terms. 



Frontiers in Cellular Neuroscience 



www.frontiersin.org 



August 2014 | Volume 8 | Article 246 | 18 



